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Abstract 

Aggregation of zero-valent nanoparticles in groundwater is influenced by several physical phenomena. The 
article shortly introduces preceding works in modeling of aggregation of small particles including influence of 
sedimentation, velocity profile of water, heat fluctuations, and surface electric charge. A brief description of 
inclusion of magnetic forces into the model of aggregation follows. Rate of influence of the magnetic forces on 
the aggregation depends on the magnitude of magnetization of the particles, radius of nanoparticles, size of the 
aggregates, and their concentration in the solution. Presented results show that the magnetic forces have 
significant influence on aggregation especially of the smallest iron particles. 



Introduction 

Zero-valent iron nanoparticles (nZVI) composed of 
iron and its oxides are spherical particles with dia- 
meter approximately 50 nm and with a large specific 
surface. These particles are used for the decontamina- 
tion of groundwater and soil and especially for the 
decontamination of organic pollutants such as haloge- 
nated hydrocarbons [1]. Nanoparticles migrate through 
the soil and can reach the contamination in-situ. Prop- 
erties of the nZVI and remediation possibilities depend 
on methods of production [2]. At the Technical Uni- 
versity of Liberec, experiments with iron nanoparticles 
TODA produced by the company Toda Kogyo Corp. 
[3] and with the nanoparticles NANOFER, produced 
by the company NANO IRON s.r.o. [4], are made. 
During a remedial intervention, transport of the iron 
nanoparticles is slowed down due to rapid aggregation 
of them. The rate of aggregation increases with grow- 
ing concentration of particles in solution and with 
growing ionic strength of the solution [5]. For the pre- 
servation of the transport properties, it is advisable to 
stabilize the particles. A lot of methods of stabilization 
were published [6-10]. We simulate the transport of 
the iron nanoparticles and that is why we examine the 
interactions among them causing the aggregation. 
Models of aggregation of small particles were pub- 
lished in many articles (e.g. [11-13]). They are mostly 
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based on the publications [14,15]. However, this gener- 
ally used model is insufficient for our case. A surface 
charge established on the surface of particles causes 
repulsive electrostatic forces between them. The influ- 
ence on the aggregation into the known aggregation 
model was implemented [16]. The iron particles cor- 
rode in the water, and this process causes change of 
the surface charge as well as the change of the rate of 
aggregation [17]. Because the particles are made from 
iron, they also have magnetic properties, which signifi- 
cantly affects the rate of aggregation [2,18-21]. That is 
why we want to derive a mathematical model of mag- 
netic forces among particles and to add it into the 
aggregation model. There is shown a procedure of the 
derivation of the model in this paper and there is 
made also an evaluation of the rate of aggregation 
influenced by magnetic forces here. 

The extended model of the aggregation of iron nano- 
particles will be included into a solver of particle trans- 
port in groundwater. It would allow to simulate the 
transport of iron nanoparticles and to predict the effi- 
ciency of the remedial intervention. That could be useful 
for the proposal of optimal remedial intervention, which 
would enable to decontaminate an affected area effi- 
ciently and economically. 

Aggregation of Colloids and Small Particles 
in Groundwater 

The particles in groundwater aggregate easily. They cre- 
ate clumps of particles up to the size of several urn [20] 
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that cohere and decrease the possibility of migration of 
particles through pores of the ground. The aggregation 
of particles is proven by experiments described in many 
articles. In [22], particle size, size distribution and sur- 
face composition were characterized by transmission 
electron microscopy (TEM), X-ray diffraction, high- 
resolution X-ray photoelectron spectroscopy, X-ray 
absorption near edge structure, and acoustic/electroa- 
coustic spectrometry. There are presented micrographs 
of a single particle and aggregates of iron particles in 
the article. In [3], characterization of iron nanoparticles 
using TEM according to methods of its preparation is 
done. In [20], a type of aggregation according to begin- 
ning concentration of iron nanoparticles is studied by 
dynamic light scattering, optical microscopy, and sedi- 
mentation measurements. 

The aggregation of the particles is caused by many 
processes that generally occur during the particle migra- 
tion. The decrease in the mobility can be formulated by 
a rate of aggregation that is given by the mass transport 
coefficients |3 [m 3 s" 1 ]. It was published in many papers 
(e.g. [11,15]). The coefficients give a probability Py of 
creation of aggregate from particle i and particle / 
together with concentrations «„ «, of particles i and par- 
ticles / (1). Particle i means the aggregate created from i 
elementary nanoparticles. 



PijUiTlj, 



(1) 



(2) 



The pfj is the mass transport coefficient of aggrega- 
tion caused by the gravity, /?| is the mass transport 
coefficient for the velocity gradients, and stands for 
the mass transport coefficient of the heat fluctuations. 
The notation is adopted from [11]. 

The first process that causes the aggregation of parti- 
cles is sedimentation. Due to the gravitation forces, the 
particles fall. According to their size, the velocity of the 
sedimentation is different for the different aggregates. It 
implies that the particles have a chance to meet the 
others and aggregate due to the attachment forces. Fol- 
lowing [11], the mass transport coefficient for sedimen- 
tation is 



72?7 



{ Qt -Q){d i +d i ) 2 \df 



(3) 



where g is gravity acceleration, r| is viscosity of med- 
ium, Q is density of medium, Q p is density of aggregating 
particles, d t is diameter of particle i. 

The process that works similarly is water drifting. The 
flowing water in a pore in soil has a velocity profile and 



in the middle of the pore the velocity of flowing is lar- 
gest. Since the particles have different velocities, accord- 
ing to the place where they are in the pore, the particles 
can meet each other and create the aggregate. Following 
[11], the mass transport coefficient for the velocity gra- 
dients of particles is 



p\=-G{d i + dry 



(4) 



where G is the average velocity gradient in a pore. 

In the case of the small particles, the heat fluctuation 
of particles has a significant effect on the particle aggre- 
gation. Brownian diffusion causes a random movement 
of the particles and again it facilitates the aggregation. 
Following [11], the mass transport coefficient for the 
Brownian diffusion is 



Pi 
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where I<b stands for Boltzmann constant and T 
denotes absolute temperature. 

A statistical assessment of the importance of the parti- 
cular processes to the creation of the aggregates was 
done. The Table 1 shows that for the smallest particles 
the Brownian diffusion is most considerable. The sedi- 
mentation is most significant when the difference 
between sizes of the aggregates is largest. It is a conse- 
quence of the fact that the difference between the velo- 
cities of the particles is largest. The velocity gradients 
depend on the pore size. When the size of pores is 
small, the aggregation is most influenced by the velocity 
gradient, if the difference between the particles is large. 

The mass transport coefficient for the velocity gradi- 
ents is quantified for the case with a small size of pores 



Table 1 The mass transport coefficients for Brownian 
diffusion, velocity gradients, and sedimentation, for 
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and a large flux (flux = 3.67 x 10" ms" , porosity = 0.39, 
velocity gradient G = 50 s" 1 ). In the other cases, the 
mass transport coefficient would be much smaller than 
the others. 



particles. The total probability of the aggregation of a 
particle i with a particle / is then 



(Pl el + Pl el + Pl el )n inj . 



(12) 



Electrostatic Forces Among Iron Nanoparticles 

A mathematical model of aggregation for the case of 
iron particles was compiled. The sedimentation, velo- 
city gradients, and Brownian diffusion are not suffi- 
cient for the description of the process aggregation. 
The surface of iron particles oxidizes. The ions on the 
surface attract ions from the electrolyte and the elec- 
tric double layer arises. Therefore, the influence of the 
electrostatic forces was added in the mass transport 
coefficients. The detailed derivation of it is published 
in [16,23]. 

For sedimentation, the mass transport coefficient /Jjj 
has the following form: 



,3d 



Pi' 
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where a ; and a, stand for surface charge on particle i 
and particle j, respectively, e 0 is permittivity of the med- 
ium. If the term that reduces the mass transport coeffi- 
cient is greater than the mass transport coefficient 
without the influence of electrostatic forces /?| , the 
probability of collision of particles i and /' should be 
equal to zero. That is why 



~Pf 



(7) 



For the velocity gradients, the mass transport coeffi- 



cient fif: el is equal to 
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And finally, the Brownian diffusion gives the mass 
transport coefficient ji 



lei 



pf=p\- 



ndfdjdjCjj 
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Plf =max(0,Af). 



(10) 



(11) 



The probability of particle collision decreases quadra- 
tically with the quantum of the surface charge of 



Values of mass transport coefficients for aggregates 
with the sizes between 50 nm and 5 um are com- 
puted. The surface charge depends on £ potential (see 
e.g. [24]). £ potential depends on pH of the water. The 
measured results of this dependence acquired using 
the Malvern ZetaSizer are shown in the Figure 1. 
Zero-valent iron provides alkaline reaction in water, 
so the measurement was done for higher pH values 
only. 

The statistical assessment of the importance of the 
electrostatic forces to the creation of the aggregates was 
done in [16,23]. The results are that the mass transport 
coefficient for Brownian diffusion is limited by the elec- 
trostatic forces mostly for large aggregates. The mass 
transport coefficient for the velocity gradients is not lim- 
ited by the surface charge for measured £ potential. The 
t, potential would have to be minimally ten times larger 
to affect the rate of aggregation. The mass transport 
coefficient for sedimentation is limited by the surface 
charge mostly for the small particles. 
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Figure 1 Dependence of ; potential of nZVI on pH measured 
with Malvern ZetaSizer 
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Magnetic Field of Iron Particles 

The iron particles NANOFER that are used for the reme- 
dial intervention were measured by magnetometer MPMS 
XL, an equipment based on the SQUID effect (Supercon- 
ducting quantum interference device), owned by the 
Palacky University Olomouc, Czech Republic. The iron 
particles are ferromagnetic, a hysteresis loop of the iron 
particles measured by SQUID is in Figure 2. That is the 
reason why it is necessary to include the influence of the 
magnetic forces among particles to the model of aggrega- 
tion of the particles. However, it is a very complicated pro- 
cess which cannot be described analytically. That is why 
only the description of the magnetic forces between two 
particles is shown in this paper. Although it is not a suffi- 
cient model, it can be used for determination of something 
what we could call "effective range" of the magnetic forces. 
Also the limits of the sizes of the magnetic forces can be 
determined. It can be used for the assessment of the aggre- 
gation of the particles depending on their concentration. 

According to [25], the electromagnetic potential in the 
point r near a permanent magnet is equal to 



CM 
J R 



MR 



V, 



(13) 



where the vector M is the vector of magnetic polariza- 
tion at the point dV, the vector R is the difference 
between the source of magnetic field dV and the point 
r, R is the length of R. 



Intensity of the magnetic field H can be subsequently 
computed as 



H(r) = -grader)). 



(14) 



Finally, the magnetic force between the source of the 
intensity of magnetic field H and a permanent magnet 
with the vector of polarization M 0 at the point r is 
equal to 



F(r) = - J(M 0 ■ grad)H(r)dV. 



(15) 



Magnetic Forces Between Two Spherical Iron 
Nanoparticles 

The scalar potential of the magnetic field around one 
homogeneous spherical iron nanoparticle with radius a 
located at the point [0, 0, 0] was determined: 



(* 3 - rcomv ™m 4/ded(Pi (16) 



2.71 K a 



where a is the radius of the nanoparticle and x 2 , 
x 3 ] are the coordinates of the point r. The direction of 
the vector of polarization M is set to the direction X3, M 
is the magnitude of the vector M. 

After integration, the magnetic potential around a fer- 
romagnetic sphere is obtained: 
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Figure 2 Hysteresis loop of the zero-valent iron nanoparticles 
for temperature of 300 K measured with magnetometer MPMS 
XL by Jiff Tucek at the Palacky University Olomouc 
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According to (14), the components of the vector of 
intensity of the magnetic field around a spherical ferro- 
magnetic particle is 
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where 5, 3 stands for Kronecker delta and i = 1, 2, 3. 
The symbol C(r) replaces 



C(r) = tJxi + x\ + x\ - a 2 . 



(19) 
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The derived formula of the size of the magnetic forces 
between two iron nanoparticles is very extensive; hence, 
it is not presented here. Though, an example of the 
numerical result is shown. In Figure 3, there is the 
visualisation of a part of the vector field of the magnetic 
forces between two nanoparticles. First nanoparticle is 
in an arbitrary point near second nanoparticle with 
radius a which is touching the center of the upper right 
side of the figure. The figure is created by the software 
Mathematica 5, copyrighted by Wolfram Research, Inc. 

Magnetic Field Around an Aggregate 

The aggregate of iron nanoparticles is in fact a clump of 
many permanent magnets. It is impossible to establish an 
analytical model of interaction of two such aggregates. 
To analyze statistically the influence of the magnetic 
forces on aggregation of two nanoparticle aggregates, a 
script for the examination of the worst possibility (the 
largest forces) and the averaged possibility of influencing 
the aggregation by the magnetic forces was written down. 

The statistical model of the aggregate is made so that 
the volume of the aggregate is filled by uniformly distrib- 
uted nanoparticles (small homogeneous magnets) with 
randomly uniformly distributed direction of the magnetic 
polarization. The magnitude of the magnetic polarization 
of all nanoparticles is the same. The magnetic potential 
of the aggregate is then the superposition of magnetic 
potentials of all nanoparticles creating the aggregate: 
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Figure 3 Visualization of the vector field of the magnetic 
forces between two spherical particles of nZVI, using software 
Mathematica 5, copyrighted by Wolfram Research, Inc One 

nanoparticle is in an arbitrary point near a nanoparticle with radius 
a which is touching the center of the upper right side of the figure. 



where fa is potential of magnetic field of the nanopar- 
ticle located in the point r t . 

The influence of magnetic forces in comparison with 
the gravitational forces is investigated. It could be com- 
pared also with other affecting forces but the gravitation 
force was chosen for the reason of small number of 
variables. If one aggregate is in a fixed position and 
another one is located somewhere vertically under it, 
there should be a unique distance ("effective range") 
between them so that if they are closer than it, the 
lower aggregate would attach to higher one by magnetic 
forces. If they are more distant, they would sediment 
separately. The distance is measured between the sur- 
faces of particles. 

In Figure 4, the dashed line characterizes the effective 
ranges for interaction of chosen aggregates interacting 
with a single nanoparticle. The solid line in the graph 
characterizes the effective range for the interaction of 
two aggregates of the same size. The absolute value of 
the magnetic force and consequently also effective range 
quadratically depends on the magnitude of the magnetic 
polarization. The graph is plotted for the magnetic 
polarization 170 emu g' . More important information 
than absolute values is the trends of the lines. 

According to the results from the Table 1 in the Sec- 
tion 2, sedimentation influences the aggregation mostly 
when the difference between sizes of the two aggregates 
is large. Consequently, the dashed line in the Figure 4 
gives a good information about the influence of mag- 
netic forces on aggregation of aggregates of different 
sizes. The solid line comparing the influence of mag- 
netic and gravitational forces between two similar aggre- 
gates does not include the real information about the 
influence of magnetic forces because another force than 
gravitational governs the aggregation process in such a 
case. 
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Figure 4 Effective range of the magnetic forces of chosen 
aggregates. The dashed line characterizes the effective ranges for 
the interaction of aggregates interacting with a single nanoparticle. 
The solid line characterizes the effective ranges for the interaction 
of two aggregates of the same size. The graph is plotted for the 
magnetic polarization 170 emu g , 
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On the basis of this result, it is obvious that the mag- 
netic forces among particles have significant effect on 
the rate of aggregation of the particles. On the basis of 
the effective range, concentration of particles uniformly 
dispersed in solution can be computed. According to 
distances between particles, the concentration would 
have to be very small (about 15 mg 1" ) to be possible to 
neglect the influence of the magnetic forces. 

Conclusion 

The influence of the magnetic forces among iron parti- 
cles on the rate of aggregation in terms of the "effective 
range" was assessed. The effective range is a distance in 
which the magnetic forces outweigh the gravitation 
force that causes the aggregation. To assess the mag- 
netic field around more interacting aggregates, it should 
be made a model using the Finite Element Method 
(FEM) or another numerical method. The next steps in 
the studying of nanoparticle aggregation due to mag- 
netic forces can be the evaluation of the effective range 
in comparison with other forces indicated in this paper 
and/or building a FEM model of nZVI aggregate. 
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